Effects of coating rate on morphology of copper surfaces 
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Abstract 

We have used standard fractal analysis and Markov approach to obtain further insights on 
roughness and multifractality of different surfaces. The effect of coating rates on generating topo- 
graphic rough surfaces in copper thin films with same thickness has been studied using atomic force 
microscopy technique (AFM). Our results show that by increasing the coating rates, correlation 
length (grain sizes) and Markov length are decreased and roughness exponent is decreased and our 
surfaces become more multifractal. Indeed, by decreasing the coating rate, the relaxation time of 
embedding the particles is increased. 
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I. INTRODUCTION 



Surface, which is the first interface of a material, has an important role in the interaction 
of matter with environment. Physical and chemical properties of the surfaces are not only 
determined by the material properties but also to a significant degree by the topography. No 
surface is perfectly flat and every surface depending on the scale it is observed, has a certain 
amount of roughness. Surface roughness is one of the important properties of the surfaces 
that affect many other properties of the surface like adhesion, friction and contact, reflection 
or scattering [lKD]. The roughness effect appears in various devices such as field emission 
devices js], sensors js|, self-cleaning materials [10]. It should be noted that roughness is not 
an intrinsic property of the surface. Indeed, it depends on the scale of observation. In other 
words, when we observe a surface from different scales different roughness could be obtained. 

In this article, we investigate the effect of coating rate on statistical properties of our 
prepared surfaces. These surfaces have the same composition and have been prepared in the 
same condition. Coating rates of these samples differ while final thickness of the samples is 
the same (250nm). We have considered copper thin films. Copper thin films have different 
applications in various technologies such as optics and laser science, because of the high 
reflection power of copper in red and infrared region of spectrum. Glass substrates were used 
in this survey because of smoothness. In Tab. I, the experimental conditions of two selected 
thin layers are given. AFM measurements of these samples were carried out and the exported 
data was used for further calculations. The topography of the samples was investigated using 
Park Scientific Instruments (model Autoprobe CP). The images were collected in a constant 
force mode and digitized into 256 x 256 pixels with scanning frequency of 0.6 Hz. A rough 
surface can be described mathematically as h(x), where h(x) is the surface height of a rough 
surface with respect to a smooth reference surface defined by a mean surface height and x 
is the position vector on the surface. 

To study the effect of topography and scaling properties of the thin films, the standard 
fractal analysis is used [ill Q]- When long-range correlations are absent in h(x), short- 
range correlations may exist. In this case, it can be more suitable to study multifractality 
by Markov analysis [bj, [ijj . Fazeli at al. showed how Markov processes play a fundamental 
role in probing rough surfaces and characterizing their topography. They explained tip 
convolution in AFM images by non-Markovian properties in the AFM images reported. Our 
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FIG. 1: (Color online) Atomic force microscopy (AFM) surface images of two samples: (a) sample 
1 represents a surface with roughness 5.2nm and (b) sample 2 represents a surface with roughness 
3.6nm. 

results show that by increasing the coating rates, correlation length (grain sizes) and Markov 
length are decreased, roughness exponent is decreased and our surfaces get more multifractal 
which means that our probability density functions get more non-Gaussian. 

The paper is organized as follows. Standard fractal method and Markov analysis are 
described in Section II. Data description and analysis based on these methods for two selected 
copper thin films are given in Section III. In Section IV, we present our conclusion. 



II. STATISTICAL QUANTITIES 



A. standard fractal analysis 

For a sample of size L, roughness is defined by w(L) = (((h — h) 2 )) 1 ^ 2 [11], where (• • ■ ) de- 
notes an spatial averaging samples, respectively. For simplicity, without losing the generality 
of the subject, we can assume that the mean height of the surface is zero, h — 0. 

Roughness is one of the scaling properties of the surface. Roughness scales by size of 
the system as w(L) ~ L a , where a is the roughness exponent. The common procedure to 
measure the roughness exponent of a rough surface is based on a second moment of height 
difference function defined as S 2 (l) = (\h(x + l) — h(x)\ 2 ). This is equivalent to the statistics 
based on the height-height correlation function C(l) =< h(x + l)h(x) > for stationary 
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surfaces, i.e. S 2 (l) = 2(u> 2 — C{1)). The second order height difference function S 2 (l), scales 
with I as where a = £2/2. 

Assuming statistical translational invariance, different moments of the height difference 
functions S q (l) =< \h(x + I) — h(x)\ q >, (qth moment of the increment of the rough surface 
height fluctuation h(x)) will depend only on the space difference of heights /, and has a 
power-law behavior if the process has the scaling property 

S q (l) =< \h(x + l) -h(x)\ q >oc S q (L ) I — J , (1) 

where L is the fixed largest length scale of the system, < • • • > denotes a statistical average 
(for non-overlapping increments of length I), q is the order of the moment (we take here 
q > 0), and £ q is the exponent of the height difference function. The main property of a 
multifractal process is its characterization by a non-linear £ g function of q. Monofractals are 
the generic result of the linear behavior. For instance, for Brownian motion (Bm) £ g = q/2, 
and for fractional Brownian motion (fBm) £ g oc q. 



B. Markov analysis 



For better understanding the multifractality features, it can be useful to investigate the 
multifractality by Markov analysis. When long-range correlations are absent in h(x), short- 
range correlations may exist. In this case, it can be more suitable to study multifractality 
by this approach. As a measure of surface roughness and to check the multifractal nature 
of rough surfaces, we check the Markovian nature of the increments which is defined by 
hi(x) = h(x + 1) — h(x) depending on the length scale I. 
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22? 1. If so, we es- 



First, we check whether hi(x) represents a Markov process [13, 
timate the Markov Length scale (ML); the minimum length interval over which hi(x) can be 
represented by a Markov process. For a Markov process, knowledge of P(h 2 , h\hi, l\) is suf- 
ficient for generating the entire statistics of hi(x), encoded in the n-point PDF that satisfies 
a master equation that, in turn, is reformulated by a Kramers-Moyal (KM) expansion, 



d °° / d \ k 

—p(h,l\h ,l ) = \gh) [ D k(h,l)p(h,l\h Q ,l )}. 
k=i ^ ' 



(2) 



In order to obtain the drift {D\) and diffusion coefficient (D 2 )) for 



a well defined way like it was already expressed by Kolmogorov 
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q. (2) we proceed in 



20|-|22|]. The conditional 



moments Dk(h, I) for finite step sizes Ax are directly estimated from the data via moments 
of the conditional probabilities. 

D k (h,l) = jim o T|^ - / (ti - h) k p(h',l - Ax\h,l)dh'. (3) 

For a general stochastic process, all the KM coefficients may be nonzero. However, pro- 
vided that D±{h, I) vanishes or is small compared to the first two coefficients [19], truncation 
of the KM expansion after the second term is meaningful in the statistical sense. For our 
samples, D±{h,l) is two orders of magnitude less than D 2 (h,l). Thus, we truncate the KM 
expansion after the second term, reducing it to a Fokker-Planck (FP) equation. According 
to the Ito calculus 19], the FP equation is equivalent to a Langevin equation, 

^h(l) = D^h, I) + v^(M)/(0, (4) 

where /(/) is a random force with zero mean and Gaussian statistics, 5 -correlated in h, i.e., 
< /(/)/(/') >= 25(1 — I'). Di is drift and D 2 is diffusion coefficients. To better clarify these 
variables (D\ and D 2 ), we can pay attention to the Langvin equation (Eq. 4). Here, D\ 
indicates an average height difference by walking on the surface (5h = D\5l) and D 2 plays 
the role of the variance (uncertainty) of this height hanging (f(l)D 2 5l = 5h). To establish a 
relation between Markov approach and fractal distributions we start with the Fokker-Planck 
equation for PDF of the height increments 

with the corresponding Langevin equation given by Eq. (T4J). It can be shown that for any 
series with the type of the correlations that is described by the self-affine distributions, the 
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drift and diffusion coefficients of the increment series are given by 

Dx(h,l) ~ -Eh 
D 2 (h,l)~bh 2 (6) 

where H is the Hurst exponent which refers to the first moment exponent and b indicates 
the strength of the mult i-f r act ality. It means if b = in this case we find the samples mono- 
fractal. Thus, using Eqs. §J5§, we obtain the evolution of the moments of height difference 
function, S q (l) = (\Deltah(l)\ q ) = (\h{x + 1)- h(x)\ q ), as follows 

-l^(\Ah(l)\") = q(\Ah(l)r l D 1 (Ah,l)) 

+q(q-l)(\Ah(l)r 2 D 2 (Ah,l)} (7) 




FIG. 2: (Color online) Height difference functions of two samples (1 (w = 5.2nm) and 2 (w = 
3.6nm)) for different moments q's ((a) q = 0.2, (b) q = 1.4, (c) q = 2 and (d) q = 2.4), versus 
distance. 

Then, by substituting Eq. (6) in Eq. (7) the scaling behavior of the moments of the 
height difference function are written as S q (l) ~ (l)^ q then, 

Z q = Hq-bq{q-l), (8) 

which establishes a direct link between the scaling exponents £ 9 and the results obtained for 
the drift and diffusion coefficients. 
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FIG. 3: (Color online) Drift and diffusion coefficients of the two samples (1 (w = 5.2nm) and 2 
(w = 3.6nm)) versus height which is normalized by roughness (w). 

III. DISCUSSION AND RESULTS 

In order to check the effect of coating rate on statistical properties of surfaces we have 
used samples of the same material (copper thin films) and equal thickness (250nm) which 
have been prepared in the same condition except their coating rates. We have chosen two of 
these samples (1 and 2) typically for drawing the plots. In table I, experimental conditions 
of these selected thin layers are given. AFM method was used for obtaining microstructural 
data from the surfaces (Figs. 2a and 2b). Scaling properties were obtained using AFM 
images. In order to better explain the scaling properties of these surfaces, two points should 




FIG. 4: (Color online) Exponents of the qth moments of the height difference function of the two 
samples (1 (w = b.2nm) and 2 (w = 3.6nm)) vs q. The errors from markovian methods are shown 
with the same color as the curves but the errors from the direct method are shown with different 
colors (black (sample 1) and orange (sample 2). 

be noted. 

(a) Roughness is not an intrinsic property of the surface. Indeed, it depends on the 
scale of observation. In other words, when we observe a surface from various scales different 
roughness could be obtained. 

(b) All information about the roughness is not included in the second moment of the 
height difference function [121 ] . They used the concept of the higher moments to obtain 
further insights on roughness. They showed how these concepts could explain a rough 
surface for the purpose of suitable application and how experimental parameters can affect 
its properties. 

Roughness of the two selected samples is different because of different coating rates. 
Normally, roughness is denned by root mean square of the height fluctuation (rms), however, 
roughness is not only determined from the second moment of the height fluctuation difference 
but also the other moments play role in the determination of the roughness which we call 
"generalized roughness". Fig. |2] (a-d) typically presents qth moment (0.2, 1.4,2 and 2.4) of 
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the samples' height fluctuation versus distance. Fig. [2] (c) presents second moment of the 
samples' height difference function. The slope of each curve yields the roughness exponent 
of the corresponding surface. The scale of saturation limit in this curve is the correlation 
length. At scales larger than the correlation length the second moment saturates to a height 
2w 2 . By increasing the order of moment, q, we see that the height of saturation of the 
moment is increased which is an indication of the generalized roughness. As it can be seen 
from the figure, after a particularhubi moment (q > 2) the height difference function of the 
first sample lies above the second one which represents that the concept of roughness that is 
obtained from higher moments are larger for the sample 1. In other words, higher moments 
focus on large height differences that occur in large scales. Thus, in large scales, sample 1 
is rougher than the second sample. It should be noted that the reported roughness in Tab. 
I presents the roughness for the largest scale which is the size of the system. 

When we have fractal systems, their £ g is a linear function of q and this means that they 
are the same in their scaling behavior but for systems which are multifractal this behavior 
changes and the behavior of £ g deviates from the simple linear relation with respect to q. 
As can be seen in Fig. |2j the scaling length, which is used to find the exponents £ 9 , is small. 
Thus, in order to find the behavior of the exponents of the height difference function moments 
we have used the Markov analysis too. In addition, this method helps us to better understand 
the multifractality of the surfaces under study. To use this approach we estimate the Markov 
Length scale (ML), the minimum length interval over which hi(x) can be represented by a 
Markov Process. To estimate the Markov length we have used ChapmanKolmogorov test. 
More detailed discussions of this test can be found in the appendix of [14]. The drift and 
diffusion coefficients, D\ and D 2 , were estimated directly from the data. They are well- 
represented by the approximates, 

( -0.51 ± 0.01ft, sample 1; 
Di(W) = < (9) 
[ -0.26 ± 0.03/i, sample 2. 

and 

f 0.13 ± 0.01/i 2 , sample 1; 
D 2 (h,l)=\ ' (10) 

I 0.10±0.01/i 2 , sample 2. 
Comparing Markov length scale and correlation length, ML is obtained directly from the 
joint PDF but the correlation length is obtained from height-height correlation function, 
which is the distance at which the correlation function falls to 1/e of its initial value. The 



Markov length scale of samples 1 and 2 are 1 pixel (7.8nm) and 2 pixel (15.6nm), respec- 
tively. Correlation length scale of samples 1 and 2 are 30nm and 47nm, respectively. We 
summarized these quantities in Tab. I. For sample 2, since its Markov length is other than 
one, we have considered the size effect for diffusion coefficient. The calculated diffusion co- 
efficient for this sample is the one with the correction term. Non-negligible corrections have 
oyed in order to get reliable estimates of diffusion coefficient for finite Markov 
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length, A |23|. 

D 2 (h,l) - (A/MM)) 2 , . 

° 2{K l > A) - 2A(l + ADi(M)) ' ( 5 

By obtaining the drift and diffusion coefficients (Fig. |3]), the £ q behavior from both 

standard fractal and Markov analysis has been calculated. Fig. H] shows the results from 

both methods. The curves with error-bars are the ones calculated from the Markov analysis 

and the other two curves are from standard method. The equation for the £ q for the two 

samples from the Markov analysis are £(g) = (0.26 ± 0.03)g — 0.10 ± 0.01q(q — 1) for sample 

1 and £(q) = (0.51 ± 0.01)g - 0.13 ± 0Mq(q - 1) for sample2. 

As can be seen in Fig. HJ £ 9 of both of our samples deviate from lines and the one which has 
larger root mean square (sample 1) deviates more. Higher moments have more information 
from the larger fluctuations, which in correlated systems larger fluctuations appear in the 
larger scales. Here, £ g of our samples shows how large scales play more effective role in the 
height fluctuation of sample 1. This means that large moments of the rougher sample grow 
more with respect to the other sample. In other words, as can be seen in Fig. [2j in frame 
(a) the q = 0.2 moment of height difference function of sample 1 is below sample 2 (after 
saturation length), but by increasing q this behavior changes. Lower moments of sample 2 
are higher than sample 1 and in the higher moments sample 1 is higher. Lower moments 
describe the small height fluctuations in the surface. 

To summarize, by increasing the coating rates, correlation length (grain sizes) and Markov 
length are decreased and roughness exponent is decreased and our surfaces become more 
multifractal. 



IV. CONCLUSION 

Different conditions in particle deposition leads to the change in the topography of the 
surfaces. We have shown that by changing the coating rate properties of the surface differ. 
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TABLE I: Experimental conditions and properties of the two selected copper thin layers: coating 
rate, coating time, thickness (r), Roughness (w), pixel's size, correlation length (l CO r), roughness 
exponent (a) and Markov length scale. 



Sample 


Coating rate (2*2) 


t(sec) 


r(nm) 


w(nm) 


Pixel's size (nm) 


l cor (nm) 


a 


ML (nm) 


1 


1.9 ± 0.04 


128 ± 1 


250 ± 5 


5.2 ± 0.2 


7.8 


31 ± 7.8 


0.20 ± 0.01 


7.8 


2 


1.7 ± 0.04 


145 ± 1 


250 ± 5 


3.6 ± 0.2 


7.8 


47 ± 7.8 


0.40 ± 0.01 


15.6 



The average energy of the depositing particles is increased by raising the coating rate. In 
this case, the particles penetrate through the surface and this affects the roughness and 
multifractality of the surface. 
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